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Abstract 

A model of a high subsonic jet with a nearby ar- 
ray of flexible, aircraft-type panels is studied numer- 
ically in two dimensions. The jet is excited by a lim- 
ited duration, spatially localized starter pulse in the 
potential core. The long time evolution of unsteady 
disturbances in the jet, the responses of the panels 
and the ensuing radiation are computed. The results 
show that the spectral response of both the jet and 
the panels is concentrated in a relatively narrow fre- 
quency band centered at a Strouhal number (based 
on jet exit velocity) of approximately 0.25 and as- 
sociated harmonics. The loading on the panels gen- 
erally increases with downstream distance. Panel 
radiation is weakest in upstream directions. Interior 
zones of silence, due to destructive interference of 
radiation from the panels, are observed. 

1 Introduction 

This paper describes the results of a numer- 
ical simulation of jet noise in the presence of four 
flexible aircraft-type panels in a panel-stringer as- 
sembly. The simulation is based on a model which 
fully couples the fluid dynamics of the jet flow to 
the panel motion and the resulting acoustic radia- 
tion. The primary objective is to determine the role 
played by near sonic jet exit conditions on installa- 
tion effects from the nearby flexible structure and 
on the response and the acoustic radiation from the 
structure. 
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The long time evolution of unsteady disturbances 
in a nearly sonic jet exiting from a converging noz- 
zle, together with panel responses and ensuing radi- 
ation from the panels are considered in this paper. 
In previous work, jet acoustics, panel response and 
radiation have been considered for lower speed jets 
exiting from straight pipes, both at rest 1 and in for- 
ward motion. 2 The present work focuses on jet evo- 
lution over much longer timescales, emphasizing in 
particular vortex shedding from the nozzle lip. 

In previous analyses, the exact sources of jet 
noise have been identified from the basic equa- 
tions of fluid dynamics. 3, 4> 5> 6 Generally, these ex- 
act sources have to be modeled in some way for 
computation. An important feature of the present 
method is the direct computation of at least some 
of the natural sources of jet noise, namely fluid dy- 
namical instability waves which develop due to in- 
stability of the jet shear layer. Experiments have 
demonstrated the existence of large scale structures 
or instability waves in jets. 7, 8> 9 These structures 
are believed to act as sources of sound, a point 
also confirmed by analytical studies 10, n ’ 12, 13 and 
computations. 9, 14, 15 

The geometry of the computational model is 
shown in figure 1. The Euler equations are solved in 
two domains; the jet domain and the radiation do- 
main separated by an array of four flexible panels. 
Panel response and radiation are also computed and 
are fully coupled to the fluid dynamics in the sense 
that at each timestep the fluid dynamics (Euler) 
computation provides the pressure difference across 
the panels, thereby allowing computation of panel 
displacement and velocity. The resulting panel ve- 
locity then serves as a boundary condition for the 
Euler computation. Thus no modeling of the panel 
excitation is needed. Other approaches to the gen- 
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Figure 1: Computational Domain 


era! problem of describing panel loadings are possi- 
ble, e.g., Dowell . 1 ' 3 Similar computations for bound- 
ary layers and for panels Excited by large amplitude 
acoustic disturbances in an ambient medium have 
also been performed . 1 ' 1 18 

The paper is organized as follows. In the next 
section there is a description of the model and a 
discussion of the numerical method and boundary 
conditions. Following that the results are presented 
first, for the case of a free jet and then for the full 
jet flow/acoustic/structural interaction. Finally, the 
results are summarized and conclusions are drawn. 

2 Numerical Method 

The computational domain is shown in figure 1. 
Unsteady pressure, density and velocity are com- 
puted in two regions. The lower domain in figure 1 
contains the jet, exiting from a contracting nozzle of 
width D, while in the upper domain, (modeling the 
aircraft interior) only disturbances due to the panel 
radiation are present. These two regions will be re- 
ferred to as the jet and radiation domains, respec- 
tively. The wall boundary is a rigid wall containing 
four adjacent flexible panels (denoted as panels 1-4 
in figure 1) with rigidly clamped boundaries. The 
panels vibrate in response to excitation by sound 
from the jet, thereby radiating acoustic disturbances 
into both domains. In the jet domain, the panel ra- 
diation level is significantly smaller than the acoustic 
disturbances generated in the jet, in contrast to the 
radiation domain where panel radiation is the sole 
source of acoustic disturbances. 

The numerical method involves coupling the com- 
putation of a nonlinear equation governing the panel 


responses (the beam equation) to an Euler compu- 
tation performed in both the jet and radiation do- 
mains. The panel vibration is fully coupled to the 
fluid dynamics in that at each time step the pres- 
sure difference across the panels, computed from 
the Euler computations., serves as a forcing term for 
the beam equation. Similarly, the displacement ob- 
tained from the beam equation is differentiated in 
time and is imposed as a boundary condition for 
the Euler computations. The numerical method has 
been described in detail . 1 Thus, the presentation 
here will be brief. 

The nonlinear beam equation is 


8 4 z 0 2 z 0 2 z dz , 
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where z represents the beam transverse deflection, 
Pi is the mass per unit volume of the beam, h is 
the beam thickness, 7 is the physical damping, Dj = 
Mh 3 / 12(1 — v 2 ) is the stiffness of the beam, M is 
the modulus of elasticity and v is the Poisson ratio 
of the beam material. The coefficient N x {z) of the 
nonlinear term represents the tension created by the 
stretching of the plal< due to bending. The pressures 
in the radiation and jet domains are p + and p ~ , 
respectively. The solution of equation ( 1 ) is obtained 
at .each time step using an implicit finite difference 
method. The panels are assumed clamped at both 
ends. 

The coupling of the beam computation to the Eu- 
ler computation occurs through the forcing term 
given on the right hand side of equation ( 1 ) . The 
pressures p + and p~ are obtained from the Euler 
computation using an explicit scheme. The displace- 
ment. at. the new time level is then obtained from 
solving equation (1) one time step. The normal ve- 
locity, v, is then obtained from differentiating z and 
employed as a. boundary condition to complete the 
update to the Euler computation. Since this pro- 
cedure is employed at each time step, the fluid and 
structural calculations are fully coupled. 

The Euler equations are solved in conservation 
form for the vector 


w = {p,pu,pv,E) T , 

where p is the density, u, v are the x and y compo- 
nents of the velocity respectively and E is the total 
energy per unit, volume, 

E = §p( u" + V) + c v pT , 

where T is the temperature and c v is the specific heat, 
per unit, volume. The pressure, p, is obtained from 
the equation of st.a.t.fe,. The Euler equations are solved 
separately in both the jet. and radiation domains. 
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In the jet domain the Euler equations are modi- 
fied to account for the jet flow. The jet exits from 
a nozzle of width D and the solution is computed 
both within and exterior to the nozzle. The Euler 
equations are modified to account for two different 
non-homogeneous forcing terms. 1 One term serves 
as a starter pulse to excite the jet. It corresponds to 
a localized source of mass injection at the location 
(x s , yj ), where yj is the location of the jet axis (3.5 D 
from the wall) and x s is approximately 1.15D. An 
alternative approach, involving time harmonic exci- 
tation of the jet is described by Mankbadi et al. 19 
The second forcing term is designed so that in the 
absence of the starter pulse the solution to the Euler 
equations would be a stationary profile correspond- 
ing to a spreading jet. The inclusion of this term sep- 
arates the computation of the disturbance, in par- 
ticular the resulting instability waves, from the com- 
putation of the mean flow (i.e. the spreading jet). 
Thus, the resulting system of equations allows for 
the simulation of instability waves and the result- 
ing sound generation, together with the curving of 
acoustic waves in the jet flow Held, without requir- 
ing the computation of the spreading jet itself. Al- 
though this is a simplified model, it captures many 
of the observed features of instability wave generated 
jet sound and permits high resolution computation 
of the coupling of jet noise with the flexible pan- 
els and the resulting radiation from the panels. In 
particular, the model allows for computation of the 
natural sources of jet noise (the instability waves) 
together with the sound radiated by these sources. 

The initial conditions are taken to be ambient 
data in the radiation domain and the mean state 
in the jet domain. The model for the mean axial 
and normal velocities is obtained from Maestrello 20 
in Cartesian coordinates. In contrast to previous 
calculations 1, 2 at lower jet speeds, the current study 
requires a mean flow model which addresses the com- 
pressibility of the flow. Thus, the mean pressure 
and density are also modeled. Inside the nozzle, 
the quasi 1-D equations of isentropic flow are em- 
ployed for pressure and temperature. Outside the 
nozzle, the jet pressure is obtained from Maestrello 
and McDaid. 21 Centerline temperature is obtained 
using isentropic formulations and off-centerline tem- 
perature variations follow the axial velocity relax- 
ation. Density follows from the state equation both 
inside and outside the nozzle. It is imperative that 
the mean flow be smooth up to and including its 
second derivative for the numerical simulation of in- 
stability flows. 22, 23 Therefore, great care was taken 
in modeling the mean flow. 


The boundary conditions are as follows (refer to 
figure 1): 

1. Bounding wall - rigid conditions are imposed 
except for the flexible panels which are treated 
as described above. 

2. Nozzle exterior - impedance boundary condi- 
tions are used on the exterior of the nozzle, 
simulating the use of an absorbing material to 
absorb waves incident on the nozzle from the 
exterior. 

3. Nozzle interior - rigid conditions are used for 
the terminal two thirds of the nozzle. The 
first third (near the inflow boundary) employs 
impedance conditions to reduce the possibility 
of reflections from the nozzle inflow. 

4. Inflow for the nozzle - characteristic conditions. 
Specifically, the Euler equations are linearized 
about the ambient state, assumed to hold far 
upstream in the nozzle, and the three incoming 
characteristics are imposed. 

p + pcu, V, Coop-p/Coo 

to be the values that they would have far up- 
stream. It has been shown that this bound- 
ary condition is valid for the lowest propagating 
mode in the nozzle. 9 

5. All other boundaries in the problem are artifi- 
cial. Non-reflecting (radiation) boundary con- 
ditions are imposed to prevent spurious reflec- 
tions from propagating into the interior. These 
boundary conditions are based on a far field ex- 
pansion of the solution. 24, 25 

A finite difference scheme is employed which is 
fourth order accurate in space and second order in 
time. The scheme is a generalization of the second 
order MacCormack scheme to allow higher order ac- 
curacy in space. 26 The scheme has been discussed in 
detail. 9, 15 

3 Results 

The results are presented in two major sections; 
(i) simulation of a jet without a nearby structure (i.e. 
free jet) and (ii) the jet with the nearby structure. In 
each case, the jet exit Mach number is Mj =0.9 and 
the surrounding fluid is at rest. The nozzle utilizes 
a converging geometry with a contraction ratio of 
2:1 over a longitudinal distance of four exit widths. 
Figure 1 shows a sketch of the computational domain 
for the case which includes the structure. 
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The grids used in this study utilize stretching 
to enhance resolution in y of the shear layer and 
resolution in x near the nozzle exit. The t/-grid 
is mapped to conform to the contracting nozzle, 
while the i’-grid remains Cartesian. This devi- 
ates slightly from orthogonality in the vicinity of 
the nozzle contraction but returns to a standard 
stretched Cartesian grid in all other regions. The 
computations have been extensively validated by 
grid refinement 1 ’ 2 and, equally important, there are 
no spurious boundary reflections visible for the time 
interval considered in this paper 100). 

3.1 Free Jet Results 

For this case, non-reflecting boundary conditions 
are imposed on the upper boundary and the jet is 
centered in y. The origin of the ^’-coordinate is at 
the nozzle exit. The domain extends 45 D down- 
stream, 30D upstream and 30D in both directions 
normal to the jet axis. The computed solution is 
symmetric around the jet axis to a close approxi- 
mation, although symmetry is not imposed in the 
computation. A 581x901 grid is employed. 
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Figure 2: Near field vorticity and pressure’ contours 

The structure of the pressure and vorticity field 
near the nozzle exit are considered first. Shown in 
figure 2 are contours of the instantaneous vorticity 
and pressure fluctuation, 

p{x, y,t) = p(x, y, t) — p(x,y), 

where p denotes the mean pressure, for a small 
region near the nozzle exit at tc^/ D = 100. The 


sources simulated by the model are vorticity dis- 
turbances shed from the nozzle lip due to the near 
discontinuity in the velocity profile at the lip. The 
shedding is nearly periodic and the vortices convect 
primarily downstream. As they propagate, pressure 
disturbances associated with the vortices travel both 
upstream and downstream, ultimately leading to far 
field sound. The effect of the flow field inside the 
nozzle on the wavelength of the upstream propagat- 
ing disturbances can be seen from the shorter wave- 
lengths due to the slower upstream wave speed in 
the contracted section. 

Figure 3 shows the time history, spectrum and 
phase plot of the fluctuating pressure, p, at a point 
inside the jet flow field (x = 10D, y — y a .vis = 0.8D). 
The PSD is normalized to 0 db based on the rnaxi- 



Figure 3: Jet flow field time history, spectrum and 
phase plot at x/D = 10 and ( y — ) / A? = 0.8 

mum of the spectra. The figure shows that the spfegv 
t.ral structure of the flow field p is composed of a 
fundamental frequency, f * , such that the St.rouhal 
number, S* = /* c TC / D ~ 0.25. The spectrum also 
shows two harmonics and one quarter subharmonict 
This fundamental Strouhal number is close to obser- 
vations of peak jet output (e.g., Zaman and Yu 2 ' ). 
The precise value of S* is certainly a property of the 
particular model used in the computations. How- 
ever, this scaling for the shedding frequency has been 
confirmed for the model by computations with dif- 
ferent jet Mach numbers in a limited range around 
0.9 and with different jet widths. Examination of p 
over an earlier time window (not shown) indicates 
a more broadband structure due to sound generated 
from the leading instability wave generated from the 
starter pulse. The predominant long time feature of 
the excited jet predicted by the model is the nearly 
periodic shedding of vorticity from the nozzle lip. 
Computations with lower speed jets (not shown) in- 
dicate that this behavior becomes much more pro- 
nounced as the jet approaches a sonic exit velocity. 

The unsteady flow field in the jet generates acous- 
tic disturbances which propagate into the far field 
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as sound. Figure 4 shows a representation of p over 
the entire computational domain at tc^/D = 100. 
Although the flow is initially excited by a starter 



Figure 4: Contours of p for free jet 

pulse, at the time for the figure this initial dis- 
turbance has propagated through and exited th| 
computational domain. The remaining disturbances 
are generated due to instabilities of the jet shear 
layer. These disturbances convect downstream and 
spread with the jefc Analogous structures have been 
previously observed and computed in axisymmetric 
models.' 1 9 The figure also shows a series of waves 
propagating upstream into the far Held together with 
a similar series of waves propagating downstream 
into the far field. The waves appear to be nearly 
periodic and are of similar spatial extent. Further- 
more, they have the same frequency and wavelength. 
These waves are due to the nearly periodic shedding 
of vorticity from the nozzle lip (see figure 2) and 
the subsequent generation of pressure disturbances 
both inside and outside the jet. It can be seen from 
figure 4 that the upstream propagating waves ap- 
pear to emanate from a point closer to the nozzle 
exit than the downstream propagating waves and 
also appear to be out of phase with the downstream 
propagating waves. This results in a region of inter- 
ference between these two wave patterns along a ray 
directed approximately 100° from the jet axis. 

There is also a wave pattern within the contract- 
ing nozzle. The fluid field within the nozzle is nearly 
constant in y, leading to a pattern of nearly plane 
waves. The boundary condition inside the nozzle is 
rigid up to * = — 10 D and then is gradually changed 
to an absorbing boundary condition to help elimi- 
nate reflections from the inflow boundary. Anima- 
tion of this data indicates that only upstream prop- 


agating waves are present in the nozzle. 

The far field pressure fluctuation is plotted in fig- 
ure 5 as a function of tc^/ D at three far field angles 
(0 = 30°, 90° and 120°) taken on a circle of radius 
20 D. At all three points selected, p oscillates with a 



Figure 5: Jet far field time history, spectra and 
phase plots at 6 = 30°, 90° and 120° from the jet 
axis 

nearly constant period with a gradually varying am- 
plitude. The period is close to that of the vorticity 
generation from the nozzle lip. Thus, the frequency 
is transmitted to the far field as sound. 

In addition to the time history, the PSD of p and 
phase plots are shown in figure 5. The PSD is nor- 
malized to 0 db based on the maximum of the spec- 
tra over the till'd© angles. The figure shows that 
the far field spectral structure is similar to the jet 
flow field spectra (figure 3). The fundamental fre- 
quency corresponds to S* ~ 0.25 and the two har- 
monics and one quarter subharmonic are exhibited 
as in the jet flow. The amplitude of the first and 
second harmonics is increased, relative to the fun- 
damental, if the amplitude of the starter pulse is 
increased, confirming that these harmonics are gen- 
erated via nonlinear interactions. Furthermore, the 
presence of low level subharmonics in the pressure 
spectra, again indicative of weak nonlinearities, is 
seen. Finally, it is seen that the spectra reach a max- 
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imum at 0 ~ 30°, consistent with experiments and 
previous computations. 2, 8 ' 9 The relative harmonic 
content of the time signal increases as 9 increases. 
Thus there is a preferred upstream directivity for 
high frequencies (i.e. the harmonics). 

3.2 Jet/Structure Interaction 

Referring again to Figure 1, a wall, assumed in- 
finite in extent and parallel to the nozzle, is located 
3.5.D above the jet. The wall is assumed rigid, except 
for four flexible, aircraft-type panels with clamped 
boundaries. Each of the panels is of length 3.5D 
and thickness 0.01D. The panels are centered at 
x/D = 0, 3.75, 7.50, and 11.25, respectively. These 
are referred to as panels 1, 2, 3 and 4. Characteristic 
parameters of the panels are typical of aluminum. 

The grid in the jet domain is comparable to that 
described above for the free jet. In order to allow 
a greater distance between the jet and the artificial 
lower boundary the t/-domain extends down 45 D. 
The panels are at the origin in the vertical direc- 
tion. The Euler equations are also solved in the ra- 
diation domain. Using the fully coupled model, the 
unsteady disturbances in the jet, which serve to ex- 
cite the panels, the responses of the panels, and the 
radiation from the panels into both the jet domain 
and the radiation domain are all computed. 



Figure 6: Contours of p in jet and radiation do- 
mains 


jet. The generation of disturbances from the noz- 
zle lip is clearly visible and occurs on a continual 
basis as time evolves. There is a series of ripples 
corresponding to acoustic disturbances which radi- 
ate into the far field. The upstream radiation again 
appears to emanate from sources close to the noz- 
zle exit, while the downstream radiation emanates 
from sources further downstream resulting in an in- 
terference pattern in the range 80° <9< 100°. The 
waves radiated downstream have higher amplitudes 
and the peak radiation occurs at around 30°. 

A series of outwardly propagating ripples, repre- 
senting acoustic radiation from the panels, can be 
seen in the radiation domain. Since these ripples 
emanate from all four panels it is not possible to 
ascribe a precise geometric origin to them. Thus, 
instead of angles, far field locations in the radiation 
domain are designated by x/D values along a line 
at y/D = 25. The wavelength is comparable to the 
wavelength of the downstream radiated sound in the 
jet, i.e. the wavelength and frequency is controlled 
by S * . The wave pattern in the radiation domain 
is more sharply defined than that observed by Mc- 
Greevy et. al. 1 for a lower Mach number jet. Those 
results concentrated on earlier times and lower speed 
jets than the present computation. The largest over- 
all level of sound radiation occurs at roughly 90° and 
the level is larger in downstream directions than in 
upstream directions. 

Features of the fully coupled flow/acoustics/struc- 
tural flow field are described next. The jet domain, 
panel responses and radiation domain are examined. 



A contour plot of the unsteady flow and acoustic 
field can be seen in figure 6. In the jet domain the 
large scale structure propagating downstream is ob- 
served, trailed by other pressure disturbances in the 


Figure 7: Jet directivity of overall sound pressure 
level, with and without nearby structure 
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3.2.1 Jet Flow Domain 

In figure 7 the far field directivity is shown for 
both the free jet and the jet in the presence of a 
nearby wall. The data is given along an arc 20 D 
from the jet exit and for angles below the centerlines^ 
Far field pressures with the wall are similar to those 
described previously for the jet without any wall. 
Over the time interval considered, the differences are 
small. 

In figure 8, the incident pressure at the center of 
the panels is shown in both the time and frequency 
domains. Phase plots are also shown. The overall 
level of p increases with the downstream distance of 
the panels. The level of p on panel 1 (partially up- 
stream of the nozzle exit ) is reduced nearly an order 
of magnitude from the loading on panel 4. Since 
panel 4 is at a lower angle to the jet axis, this is in 
agreement with the directivity results in figure 7. 
The frequency spectra of the incident pressure is 


dominated by the fundamental frequency, S* j: and 
its harmonics. The relative harmonic content of the 
incident p for panel 1 is significantly increased as 
compared to the other (more downstream) panels, 
due to the preferred upstream beaming of high fre- 
quencies. 

3.2.2 Panel Responses 

The vibrating velocity, v, at the center of each 
panel is plotted in figure 9, together with the cor- 
responding PSDs and phase plots. The panels were 
held fixed until the leading pulse had passed. This 
eliminates any effect of the pulse in the panel re- 
sponses, but does lead to an impulsive start to the 
panel motion. For this reason, only the long-time 
responses of the panels are investigated. The panel 
responses increase with distance until they reach a 
maximum on panel 3 and then decrease slowly with 
distance.; The relative subharnronic level is increased 





Figure 8: Panel pressure loading time histories, Figure 9: Panel velocity time histories, spectra and 
spectra and phase plots phase plots 
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as compared to the incident pressure. This may b<|. 
due to the receptivity of the panels to low frequency 
excitation. The predominant responses of the pan- 
els, however, are at the jet fundamental frequency, 
Sf :> and its harmonics. The relatively large harmonic 
contribution to the response of panel 1 is indicative 
of the preferred upstream beaming of high frequen- 
cies, although the overall response is smaller than 
that of the more downstream panels. 

The observation made by McGreevy .ft al. 1, 2 that 
the panels act as narrow filters, converting a broad 
band input into a narrow band response, may no 
longer be relevant for jets dominated by the long- 
time vortex shedding considered here. In this case 
the panels respond clearly at the frequencies of the 
loading, somewhat enhancing the subharmonic re- 
sponse. 




Figure 10: Radiation domain far field time histo- 
ries, spectra and phase plots at x/D = —5,10,25, 
and 40 along the line y/D = '2h 


3.2.3 Radiation Flow Domain 

The panel radiation is considered next. The pres- 
sure, p, is computed along the line y = 25 D in the 
radiation domain. Shown in figure 10 are time his- 
tories, PSDs and phase plots for four different x- 
locations along this line. It can be seen from the fig- 
ure that the overall level of radiation is reduced as 
x decreases, i.e. in the upstream direction in terms 
of the jet flow field. The PSDs indicate that for 
downstream locations, there is peak radiation at 
S* , the first, two harmonics and the subharmonic. 
For upstream locations, a low level, relatively broad 
band spectrum is found. Thus the radiated sound is 
beamed primarily in vertical and downstream direc- 
tions and is controlled by S* , i.e. the vortex shedding 
frequency of the jet. 

In figure 11 the overall intensity is computed as 
a function of x. This figure illustrates the strong 



Figure 11: Radiation domain directivity of overall 
sound pressure level along the line t/ = 25D. 

vertical and downstream beaming of the panel radi- 
ation. R. also illustrates a zone of relative silence at 
about. x/D = 15 where the intensity is reduced by 
about. 15 db from nearby locations. This zone can 
also be seen in the flow field in figure 6 and is due to 
destructive interference between radiation from the 
various panels. No such zone was observed at. lower 
jet. Mach numbers 1, 2 , where the resulting radiation 
was more heavily influenced by the starter pulse and 
the sound generated by the leading instability wave 
in the jet.. 

In a. real jet., instabilities can be generated by 
bursts within the jet., generating large scale struc- 
tures, as simulated by the starter pulse and the in- 
termediate time behavior of the jet.. Also, instabili- 
ties can be generated in a. more continual fashion via. 
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vortex shedding, as is done here, by considering the 
long time evolution of the model. The model allows 
for a separation of these effects so that they can be 
studied separately. 

4 Summary 

A numerical simulation has been conducted 
of a nearly transonic jet exhausting into a fluid at 
rest both with and without the presence of a nearby 
structure. The structure is rigid except for four flex- 
ible aircraft-type panels. These panels respond to 
acoustic forcing from the jet and radiate sound into 
both the jet domain and a domain on the other side 
of the panels. The simulation utilized a full coupling 
between the Euler equations (solved on either side 
of the panels) and the nonlinear beam equation, re- 
sulting in a full coupling between the flow, acoustics 
and panel response and radiation. The results are in 
qualitative agreement with experimentally observed 
phenomena and with previous computations, indi- 
cating the present model is capable of capturing the 
important features of jet flow/acoustic/structural in- 
teractions at nearly sonic jet exit velocities. 

In particular, the results indicate the following 
properties of the long time response of excited, high 
subsonic Mach number jets: 

1. There is a nearly periodic shedding of vorticity 
from the nozzle lip. This shedding is associ- 
ated with pressure disturbances which lead to 
far Held sound. 

2. The shedding frequency scales with Strouhal 
number. S* is near 0.25 which is near the peak 
of many observed far field radiation patterns. 

3. The relative harmonic content of p increases 
with upstream angle. 

4. The far field sound peaks near 30°. 

All of the above are consistent with the features 
of previously observed real jets. Furthermore, the 
panel response and radiation results indicate that: 

1. Loading and response of the panels generally 
increases with downstream distance. 

2. The frequency spectrum of the panel responses 
is dominated by the jet shedding frequency, S* , 
its harmonics and subharmonics. 

3. The filtering action of the panels previously 
reported 1, 2 is less evident in the present work, 
where the panel forcing is dominated more by 
the periodic shedding from the nozzle lip and 
is already of a narrow band nature. 


4. The relative level of harmonics in the panel 
loading and response is increased for upstream 
panels. 

5. Radiation is primarily in vertical and down- 
stream directions. 

6. The upstream radiation is small compared to 
downstream radiation and does not appear to 
reflect the jet shedding frequency. 

7. Zones of relative silence occur due to destruc- 
tive interference. 

Finally, the code developed is computationally ef- 
ficient, running at over 210 MFLOPS on a Cray- 
YMP and over 550 MFLOPS on a Cray C-90. 
This produces a complete solution of the fully cou- 
pled flow/acoustic/structural interactions in approx- 
imately 6 CPU hours (Cray-YMP) on the relatively 
fine grids required to resolve these interaction phe- 
nomena. 
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